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I. INTRODUCTION 

Mesoscopic disordered conductors have attracted many 
theoretical and experimental investigations in recent 
years. The transport properties of chiral edge states on 
the boundary of a two-dimensional quantum Hall state 
have been of particular interest. Recently, attention 0-0] 
has focussed on a layered i/iree-dimensional (3D) system 
consisting a large number of such quantum Hall systems 
stacked upon each other. Chalker and Dohmen jl| ar- 
gued that the localized electronic states in the bulk of a 
quantum Hall system will remain localized even in the 
presence of a coupling in the third dimension, provided 
this coupling is not too strong. Electronic transport in 
such a system is therefore controlled by the itinerant, 
chiral, edge states, which now form a two-dimensional 
surface sheath on the boundary of the 3D sample. A re- 
cent experiment Q on a GaAs/Alo.iGao.gAs multilayer 
structure has indeed found that the conductivity per- 
pendicular to the layers scales with the perimeter of the 
sample when the bulk is in a layered quantum Hall state, 
thus demonstrating the existence of a conducting surface 
sheath. 

Theoretical analyses of transport in this chiral surface 
sheath have so far been carried out in the context of a 
directed network model (DN) introduced by Chalker and 
Dohmen 0. This model ignores electron-electron inter- 
actions and describes the motion of independent electrons 
along the links of a directed network; scattering events 
occur at the nodes of the network and we are ultimately 
interested in the probability distribution of transport co- 
effecients after the random scattering matrix element has 
been averaged over. The global phase diagram of the DN 
model is summarized in Fig [l]. 

We are considering a surface sheath of length L per- 
pendicular to the layered quantum Hall states, and cir- 
cumference C along the chiral edges. The system there- 
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FIG. 1. Sketch of the crossovers as a function of the length 
L. The first crossover Q is at the scale Lo = {C p a / Mq) 1 ^ 2 , 
where p s and Mq are related to microscopic length and energy 
scales in a manner discussed below (|l]) and in Ref. jjj]. The 
second crossover is near the localization length f = 8C'p s . The 
conditions for the validity of the continuum theory require 
^ S> Lo, and so the two crossover scales are well separated. 
The 2D chiral metal regime was called 0D in Ref. but we 
prefer the present terminology for reasons discussed in Ref. Q . 



fore has the geometry of a cylinder of length L and cir- 
cumference C, and we discuss the results of two-terminal 
transport measurements with leads placed at the ends of 
the cylinder. As a function of L, there are three distinct 
coherent transport regimes, separated by the lengths Lq 
and £ at which there are smooth crossovers. The precise 
values of these lengths will be discussed later, but here 
we simply note that Lo ~ (Ca^) 1 / 2 and £ ~ C, where 
a x is a microscopic length of order the spacing between 
the layers. In the presence of electron-electron interac- 
tions there will also be a temperature-dependent phase 
coherence length L v which arises from inelastic scattering 
events between electrons. This length must be accounted 
for in any comparison with experiments. In this paper we 
will, for simplicity, not consider electron-electron interac- 
tions and, therefore, our results can only be applied when 
the temperature is sufficiently low so that L v 3> L, C. 

Let us now review the basic physical properties of the 
three regimes in Fig [l]: 

(i) 2D Chiral Metal: In the plane of the layers, the elec- 
trons move ballistically in a single direction along the 
surface (the chiral edge states), while transverse to the 
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layers the motion is diffusive. As the transverse length 
L < C 1 / 2 , the electrons do not have enough time to exe- 
cute one chiral orbit in the plane before exiting through 
one of the ends of the system. The nomenclature "2D 
chiral metal" was suggested in Ref |7), from an inter- 
pretation of the structure of conductance fluctuations. 
Earlier work pj has referred to this regime as OD. 

(ii) ID Metal: The motion is as in the 2D chiral metal, 
but the system is now long enough that the electrons exe- 
cute many chiral transverse orbits before leaving through 
one of the ends of the cylinder. The probability distribu- 
tion of the transverse conductance is now very similar to 
that of a short metallic wire jjj . 

(iii) ID Insulator: If the cylinder is long enough, eventu- 
ally the one-dimensional diffusive motion along the axis 
of the cylinder undergoes Anderson localization [EJ, and 
the system behaves like an insulating wire 0. 

We now outline the results obtained in earlier studies of 
the DN model. Ref jj] presented numerical studies which 
obtained the crossover to the ID insulator. A mapping 
of the DN model to a ferromagnetic spin chain model 
was obtained in Ref ||. Ref studied the energy level 
fluctuations exclusively in the 2D chiral metal. In our 
previous work (hereafter referred to as I) , we obtained 
explicit results for the crossovers in the conductance and 
its variance between the ID metal and the ID insulator 
regimes: these results were in good agreement with the 
earlier numerical work Q. Refs Q and j(| presented an 
incomplete description of the conductance fluctuations in 
the 2D chiral metal regime, and argued these were non- 
universal. 

The present paper contains a detailed and complete 
description of the crossover between the 2D chiral metal 
and ID metal, along with explicit and exact expressions 
for universal crossover functions describing the behavior 
of the conductance and its variance. In contrast to earlier 
work HH we nn d that the variance of the conductance, 
and all of its higher moments, are universal when suitably 
expressed in terms of a small number of renormalized 
parameters which characterize the microscopic theory: 
the relationship between our and earlier results will 
be discussed in greater detail in the body of the paper. 

II. DEFINITION OF THE MODEL AND 
OUTLINE OF RESULTS 

The formalism and results of this paper are a continu- 
ation of those presented in I. We will therefore not repeat 
here the detailed discussion of the structure of the DN 
model presented in I. It was shown in I (and in Refs ) 
that the DN model can be mapped onto a continuum 
quantum field theory for a ID quantum ferromagnet in 
"imaginary time". Here the co-ordinate along the cir- 
cumference of the cylinder (of length C) is interpreted 
as the "imaginary time" direction of a one-dimensional 



quantum system whose degrees of freedom reside along 
the "space" direction of length L. The mapping involves 
an average over disorder, which is performed using ei- 
ther replica or supersymmetry methods. In the first case 
the quantum ferromagnet contains spins of the su(n,n) 
or su(2n) algebra, depending on whether we use bosonic 
or fermionic replicas, and we are interested in the limit 
n —> 0. In the second case the spins are generators of a 
superalgebra. In the present paper, we will be interested 
in the perturbative spin- wave expansion, and for this pur- 
pose it is more convenient to use a bosonic replica for- 
malism rather than supersymmetry. Although replicas 
are known to fail in some non-perturbative situations (^] 
(but not in all, see for example Ref. [0), for perturbative 
calculations they are completely equivalent to the super- 
symmetry method. If one uses, then, the bosonic replica 
formalism to treat the disorder averages, the ferromagnet 
action is (the reader is urged to consult our companion 
paper I for further details) 

Scont = _ J ^ x J dl ~ J dutrQ(u)d u Q(u)dTQ(u) 
+ y tr(^Q) 2 ) , (1) 

where Mq and p s are the magnetization per unit 
length, and the spin stiffness in the ground state re- 
spectively. Q(x, t) is a 2n x 2n matrix obeying = 
AQA, where A is a diagonal matrix with elements 
(1, . . . , 1, — 1, . . . , —1), Q 2 — I, and having n eigenvalues 
1, n eigenvalues —1, which thus parametrizes the coset 
space SU(n,n)/S(U(n)xU(n)). In the Berry phase term, 
the first term in (pj), Q(u) = Q(x, r, u) is some smooth ho- 
motopy between Q(x. r, 0) = A and Q(x, r, 1) = Q(x, r). 
The field Q satisfies the periodic boundary condition in 
the r-direction, and the following boundary conditions 
at the spatial boundary, where the DN is connected to 
ideal leads: 

Q(0,t,u)=Q(L,t,u)=A. (2) 

This paper shall focus exclusively on the properties of 
•Scout) which the reader can also consider as a ferromag- 
net of super/replica spins of interest in its own right. It 
was argued in I that the properties of the DN model are 
completely, and universally, characterized by the dimen- 
sionful parameters that appear in S^ont^ these are p s , Mq, 
L and C . Of these, L and C are given by the macroscopic 
dimensions of the sample, and therefore easily measured. 
The values of p s and M are determined by the detailed 
microscopic properties of the electronic system; neverthe- 
less if the microscopic Hamiltonian is precisely known, 
the exact values of p s and Mq can usually be determined 
exactly Q| . This is in contrast to the conventional case in 
the theory of critical phenomena, where the relationship 
of renormalized parameters to the underlying microscop- 
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ics requires solution of an intractable, strongly-coupled 
problem. 

We comment briefly on aspects of the relationship be- 
tween the values of p s and Mq for the DN model as intro- 
duced in Ref. 0, and studied in I. For the case in which 
the interior of the cylinder is supposed to be a quantized 
Hall state with Hall conductance per layer equal to 1 
(in units of e 2 /h), the parameter Mq takes the value |t]] 
l/2a x , where a x is the separation of the layers, while p s 
involves also other parameters that need not be specified 
here j7|. Here we wish to argue that, for the more general 
situation in which the bulk Hall conductance is an integer 
v > 0, say, and there are v chiral edge channels on the 
edge of each layer, the appropriate model for the large 
scale properties is again the model with action Eqn. @), 
but with Mq = v/2a x , where a x is still the separation of 
the layers. This ensures that the edge channels have the 
correct total conductance for transport along the edges, 
but the details of the distinct channels are unimportant, 
assuming that there is hopping of similar strength be- 
tween all nearby channels. If the hopping between the 
different edge channels in a layer is much weaker than 
that between corresponding channels in different layers, 
then the theoretical mappings of Refs. [|||] and I lead to v 
ferromagnetic superspin chains with weak ferromagnetic 
coupling between them. Such couplings are relevant and 
cause the behavior to cross over on large length scales to 
that for a single chain with Mq = vj2a x . (The relation 
among the values of p s for these different regimes is a 
little more complex; it can be calculated by considering 
the excited states of the coupled chains with a single spin 
flipped, in the long wavelength limit.) It is possible that 
in practise a given sample with v > 1 might not be large 
enough to reach this regime, and the comparison of ex- 
perimental data with the following theory might then be 
complicated by crossover effects. These crossover effects 
will not be further addressed here, but could be studied 
by similar methods to those below. 

There are three distinct regimes in the theory described 
by S'cont, identified by Balents et al. in Ref. M, which 
were shown earlier in Fig [j]. We can now quote the precise 
values of the scales Lq, £, at which the crossovers take 
place: 

Lq = {Cps/Mq) 1 ' 2 (3) 

and 

£ = 8Cp s . (4) 

The physics of both ID regimes was discussed in detail 
in I, where it was shown that for L > Lo, the action 
for the quantum continuum ferromagnet may be further 
reduced to that of a ID non-linear sigma model. This 
latter model has been extensively studied in the context 
of quasi-lD wires, see, in particular, Ref. p"H ]; using these 



earlier results, explicit formulas for the mean and vari- 
ance of the conductance of S'cont were given in I, valid 
thoughout the crossover between the ID metallic and ID 
localized regimes. 

In this paper we concentrate on the conduction prop- 
erties in the 2D chiral metal regime L Lq and the 
crossover to the ID metallic regime Lq <C L <C £• The 
model with action S'cont represents a continuum su(n,n) 
ferromagnet at finite effective temperature 1/C (this ef- 
fective temperature, which is really the inverse circum- 
ference, 1/C, plays a very different role than the true 
temperature, which is essentially zero, would). It was ar- 
gued in Ref |I| and I, that when the effective temperature 
1/C becomes less than the order of the low-lying level 
splittings of spin waves or magnons, which are of order 
P s /MqL 2 , there are very few thermally excited magnons, 
and the problem can be treated perturbatively. This de- 
fines the 2D chiral metal regime, 1/C <C p s /MqL 2 , or 
L -C Lq. In the ID metallic regime, the temperature 
1/C is larger than the splittings, but the thermally ex- 
cited magnons can be viewed as a slow variation of the 
spin direction with position and imaginary time. In this 
regime, the classical statistical mechanics approximation 
of neglecting the time dependence can be used, and leads 
to the ID non-linear sigma model, with coupling con- 
stant ~ 1/Cp s . Perturbation theory breaks down on 
long length scales, but is still valid when the length L 
is less than the localization length £ = 8Cp s . Thus, spin- 
wave perturbation theory can be used all the way across 
the 2D to ID crossover. In particular, the scaling forms 
of the mean and variance of the conductance can be ex- 
panded as a perturbation series in powers of L/£, times 
a universal function of L/ Lq in each term. In this paper 
we will show explicitly that for (g), the mean of g (here 
and below, the single angular brackets (...) denote the 
average of a quantity over the disorder), this takes the 
form (here and henceforth, all conductances are quoted 
in units of e 2 /h): 

(g) = £/2L + (L/0* (L/Lq) + O ((i/£) 2 ) • (5) 

To leading order in L/£ there is no dependence on L/ Lq 
(Ref. We will show below that the possible term 

of order (£/£)° vanishes identically, consistent with the 
known result that the leading "weak localization" cor- 
rection to (g) vanishes in the quasi-lD metallic regime 
L/ Lq — » oo in the present (unitary) case. The next term 
in the expansion in L/£ in (^) does have a non-trivial 
crossover at the scale Lq, described by a universal func- 
tion $ given below in Eq. (^). In contrast, for the vari- 
ance of g we obtain 

var g = (g 2 ) - (gf = $ {L/Lq) + O (L/$ , (6) 

where the universal crossover from the 2D to the ID 
metallic regime is evident in the leading term $ given 
below in Eq. (^). This leading term, which is all 
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that survives in the limit L/£ — > 0, corresponds to 
what are known as "universal conductance fluctuations" 
(UCF), which are similarly the weak coupling limit of the 
crossover function. Note that we use the term "univer- 
sal" as it is understood in critical phenomena, to mean 
that the results are independent of microscopic details 
of the model, and that the conductance fluctuations are 
described by a universal function of the sample geome- 
try (that is, in the present case, of L/Lq). This should 
be contrasted with the usage often implicit in the meso- 
scopic physics literature, in which "universal" stands for 
the value e 2 /h. In such usage, anything close to this value 
is called universal. For us, on the other hand, this unit of 
conductance is left implicit in our formulas, and it is the 
numerical coefficient, which is actually a precise function 
of L/Lq, that is universal in our sense, and need not even 
be close to unity, as in the results below. These consider- 
ations apply unchanged to the conventional UCF, which 
in our opinion are best regarded in this sense also. 

The universal functions <£> and $ are calculated in the 
following sections, and the results are 
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FIG. 2. Plots of the universal functions <E>, for the correc- 
tion to the mean conductance (top), and $, for the variance of 
the conductance (bottom). The dashed and dot-dashed lines 
are plots of the asymptotic expressions at small and large 
values of L/Lq, Eqs. ( pj]Jl2[ ) and (p|px|), respectively. 
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In these expressions A = 2tt 2 L 2 i /L 2 , and f{x) = l/(e x — 1) 
is the Bose function. These expressions are easily eval- 
uated numerically, and the resulting plots are shown in 
Fig.| 

Using the limiting behaviour of f{x) for small and large 
values of x we obtain from Eqs. (0j|) in the ID limit 
L/Lq — > oo (but still with L/£ <C 1) of long cylinders 
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These limits (p|,|To|) for long cylinders are the well known 
quasi-lD metallic results |Tl]] . In the opposite 2D chiral 
limit L/Lq — > of short cylinders 
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As an alternative to the interpretation above in terms 
of thermally excited magnons, we can also understand 
the results for this crossover in terms of the sum over 
paths taken by the propagating electrons. To calculate 
the conductance, retarded and advanced paths are re- 
quired, and these are paired up by the average over the 
disorder. This was the basis for the treatment in I. For 
(g), the leading term comes from paired paths that go 
from one end to the other of the system while propagat- 
ing a distance of order MqL 2 / p s around the circumfer- 
ence. This distance is C in the 2D regime. Therefore 
the corrections (which correspond to thermally-excited 
magnons) in the series in L/t; for (g), which reflect the 
effect of paired paths that wind around the finite circum- 
ference and interfere with themselves, are exponentially 
small in this limit. In contrast, for var g, there are non- 
trivial effects at leading (zeroth) order in which, like 
the leading term in (g), do not require winding paths. 
The correction terms to the behavior in Eq. (|l^) are also 
exponentially small in the 2D chiral metal regime. Note 
that nontrivial effects in the absence of winding paths 
also appeared in Refs. 

A result of the form of Eq. ( |l2|) for var g, has been 
claimed by Mathur J5|] and Yu [|6|, who find var g oc 
Cp s /MqL 2 in the limit L ^ Lq. Mathur's explanation of 
this result, which is essentially that the sample can be di- 
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vided into blocks of size of order L x MqL 2 / p s which have 
statistically-independent fluctuations with variance of or- 
der 1 that add to give varg, appears correct. However, 
both he and Yu || state that this result is nonuniversal, 
and Yu bases this statement on its dependence on p s /Mo 
(in our notation). Mathur's results are incomplete, as 
stated by him. His result for L <C Lo (but not that for 
L > Lo) does agree with ours in this limit for fortuitous 
reasons which will be discussed in Sec. V. Yu || states 
that his results are nonperturbative, but fail for some pa- 
rameter values. In contrast, we are able not only to com- 
pute the mean and variance of the conductance in the 
limiting cases, but also to determine the full crossover 
functions. We have already commented above on the 
meaning of universality and the scaling with L/Lq, and 
we argue that the crossover functions are universal in the 
domain of applicability of the continuum description of 
the DN model. Our calculation shows that this crossover 
from 2D to ID metal can be treated perturbatively in 
close analogy to the usual UCF in an isotropic diffusive 
system, where one can also treat the crossover from 2D 
to ID. (Here we would like to repeat a remark already 
made in I, that in the latter case, for a sample of width 
W and length L, one finds varg oc W/L as W/L — > oo. 
This behavior can be understood by viewing the sample 
as composed of independent blocks of size L x L. It is 
nonetheless universal in the sense we prefer, as discussed 
above. It is true that here there is no dependence on 
the conductivity, while there is in the chiral metal, but 
this does not mean the result for the latter is nonuni- 
versal, because the scaling functions for the chiral metal 
exhibit no-scale-factor universality, as argued in I, so the 
dependence on the bare coupling constants p s and Mq is 
meaningful.) In addition, our results differ quantitatively 
from those of Yu pi. Thus our results differ from those 
in either of Refs. [|||. 

The rest of the paper presents the details in the com- 
putation of the central results (0) and (|[). Beginning 
in the next section, we develop the spin wave expansion 
mentioned above. The su(n,n) spin chain is described in 
Sec. |n|. In the same section we introduce a parametriza- 
tion for the spins in terms of Dyson-Maleev (DM) bosons. 
In Sec. ^ we exploit gauge invariance to obtain expres- 
sions for the currents and for the moments of the conduc- 
tance in terms of DM bosons. Details of the perturbative 
calculations are given in Sec. 0. Finally, we conclude in 
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III. SPIN WAVE EXPANSION IN TERMS OF 
DYSON-MALEEV BOSONS 

In this section we will begin by considering a lattice 
discretization of iS C ont and then set up a spin-wave per- 
turbation theory using the Dyson-Maleev method. 



In I we used supersymmetry to treat the disorder av- 
erages in the mapping of the DN to a spin chain. When 
the mapping is performed instead (in a completely par- 
allel way) using bosonic replicas we obtain the following 
Hamiltonian for the spin chain: 



H 



N-l 



J(t)J(<+l) + iA(J(l) + J(JV))) . (J.-!) 



Here the J are the su(n,n) spins in a particular irre- 
ducible representation. This representation naturally ap- 
pears in the mapping as follows. We introduce n "re- 
tarded" and n "advanced" bosons a,, 6j, i = 1, 2, . . . , n. 
Their bilinear products obeying u(n,n) commutation re- 
lations are arranged in a 2n x 2n matrix 



J = 



- 2 



hSik + a\a,k 



44 



(14) 



The average over randomness in the DN model produces 



the local constraint a\a i = b[b i (we assume summation 
over repeated indices everywhere, unless stated other- 
wise) . The subspace of the Fock space of bosons speci- 
fied by this constraint forms a highest weight irreducible 
representation of the algebra su(n,n) with the vacuum 
|0) being the highest weight state. 

The trace in Eq. © is over the matrix indices of the 
spins which are multiplied as matrices. The last term 
in Eq. (|l3|) comes from the boundary, where the DN is 
connected to the ideal leads. A is the same as in (pi). 

The mean conductance (g) is given by the thermal cor- 
relator of conserved currents 



(g) = -Tr {TI(i)I{j)e 



-cn\ 



(15) 



for i 5^ j, where Tr stands for the quantum mechanical 
trace in the Hilbert space, the role of the inverse temper- 
ature is played by the circumference C of the cylinder, 
the currents are related to the spins (see I), and T 
is the time-ordering operator. Similar expressions can be 
given for the moments of the conductance. 

For the perturbative spin wave expansion it is conve- 
nient to make use of another construction of the nec- 
essary representation of su(n,n). It is very similar to 
the Dyson-Maleev construction |l^] of representations of 
su(2) and proceeds as follows. We introduce n 2 bosons, 
which we arrange in an n x n matrix b with elements 
bij, i,j = l,..-,n. To avoid confusion, we will denote 
boson creation operators by an asterisk: &*•, and reserve 
the dagger for the hermitean conjugate of matrices of 
operators. Then the hermitean conjugate matrix b' has 
elements (b' )ij — b*^. We also assume that whenever we 
write a product of several b s and b''s without indices, it 
means the usual matrix product. These n 2 bosons satisfy 
canonical commutation relations [6^,^] = 5ik$ji- Then 
the su(ro,n) spin J is given in block-matrix form by 
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J 



(n-1- btf)b n-\-btf 



(16) 



When we substitute this expression for the spins J in 
terms of DM bosons into the Hamiltonian ( |l3| ) and bring 
them to normal order, we obtain, in particular, quartic 
terms which do not have the matrix product index struc- 
ture. Next we introduce a bosonic functional integral 
using bosonic coherent states \b) satisfying b\b) = b\b). 
In this functional integral the bosonic operators b are re- 
placed by commuting variables b (without a hat), and 
we can bring the quartic terms back to the usual matrix 
product form. This gives the following action: 



S= f drVtr(6t(j)5^(j) 
Jo 



+ K (ftt(i+l) - &t(t)) (b(i + l) - b(i)) 
-Kb\i) (b(i+l) - &(»)) &t(*+i) (b(i + l) - &(»))) 

(we dropped a constant n/2 which disappears in the 
replica limit). The measure for the functional integral 
over b, w is the usual one, D[b, b']. 

Next we take continuum limit in the spatial direction 
and obtain the action S = So + S- m t, where 

S = dx rfrtr {2M Q tfdJ> + Ap s d x tfd x b) , (17) 
Jo Jo 

= -Ap a / dx / drtr (tfd x btfd x b) , (18) 
Jo Jo 

with Mo = l/2a x and p s = Ka x . The r-derivative term 
in So does not have the canonical form. This could be 
mended by a rcscaling of the bosonic fields. However, we 
will keep the present normalization to avoid cumbersome 
coefficients in different expressions appearing later. The 
boundary terms in the Hamiltonian (|l3|) in the contin- 
uum description force the field b(x) to take zero values 
at the boundaries: 



6(0, t) = b(L,r) = 0. 



(19) 



Before giving the expressions for the currents in terms 
of DM bosons and evaluating their correlators in the the- 
ory with the action ( |l7j , |l8| ) , we note that the latter could 
have been obtained from S con t of Eq. (|l|) by the following 
formal change of variables in the functional integral over 
the field Q: 



Q = 



1 + 2btf 
-2(1 + 66 + )6 



26+ 
-1 - 2btf 



(20) 



Even though this expression does not satisfy the original 
conditions on Q (stated after Eq. 1), it formally coincides 
with 2 (6 1 J 1 6) where J is given by ( |l6| ) (n in the lower 
left corner of ( ]lr^ ) disappears after normal ordering the 
term bb'). Also, the SU(n,n)-invariant functional mea- 
sure for Q reduces to D[b, tf]. These observations give 



us a very general and straightforward way of obtaining 
the expressions for the currents and the moments of the 
conductance. Namely, we have to make the action (|l|) 
gauge invariant by coupling to a gauge field. Then func- 
tional derivatives with respect to this field will give the 
currents and the moments of the conductance. When 
used together with Eq. (p0|), this procedure gives all the 
quantities in terms of DM bosons. In the next section we 
describe this procedure in detail. 



IV. GAUGE INVARIANCE, CURRENT 
CONSERVATION AND MOMENTS OF 
CONDUCTANCE 

It is clear from the last section that the distribution of 
the conductance of the DN model is related to correlators 
of the conserved current of the ferromagnetic spin chain. 
The computation of the required correlators is greatly 
simplified by an understanding of the gauge-invariance 
properties of S cont which will be discussed in this section. 

We can make the action (|l[) invariant under the local 
gauge transformation Q — > Q = UQU^ 1 , — > — 
UAfJJ- 1 - iUdfJJ- 1 , with U G SU(n,n), by replacing 
the partial derivatives with the covariant ones: d^X — > 
D^X = dfjX + i[A^,X], where the gauge potential A^ 
is an element of su(n,n). However, as discussed in detail 
in Ref. [[uj, this is only necessary for the (d x Q) 2 term. 
The Berry phase term, being a total derivative, can be 
made gauge invariant by adding a boundary term, and 
the gauge invariant action is 

S[A] = - f dx f dr f dutYQ{u)d u Q{u)d 7 Q(u) 
Jo Jo V 2 Jo 

+ % M trA T Q + P±tv{D x Q) 2 ^j . (21) 

Gauge invariance of this action implies some conserva- 
tion laws, or Ward identities (WI's). In their derivation 
we closely follow Ref. Let us introduce the generat- 
ing functional 



(22) 



Z[A] = Jd[Q] exp(-5[A]). 
Using the notation 

((X))=Z- 1 [A] [D[Q]X[Q]exp(-S[A]), (23) 



where X[Q] represents any functional of Q, we can show 
that gauge invariance of the action ( f2l| ) results in the 
following equation of motion: 



((SS/SQ)) = 0. 



(24) 



As shown in Ref. EJ, this is equivalent to the following 
covariant current conservation law: 



G 



Dr((jr))+D x (( ]x ))=Q., 



(25) 



where we introduced matrix-valued gauge invariant cur- 
rents with elements 



(55 [A] 



(26) 



Equation ( p5| ) is the first in the series of WI's coming 
from the gauge invariance. The gauge invariant currents 
(^6|) contain gauge potentials, and further WI's are 
obtained differentiating with respect to A. In particular, 
we can obtain such identities for the mean two-probe 
conductance and its variance, which in this field-theoretic 
formalism are obtained as follows. Let us assume that the 
source field A is independent of t, and moreover, has the 
special form 



A T = 0, A x (x) = 



A«(a:) A+(x) 



A-{x) A^\x) 
Then the mean conductance is given by 

S 2 Z[A] 



(27) 



(g(x,x')) 



lim 



SAt 1 (x)SA u {x f ) 
If we introduce the spatial currents integrated over r, 
6S[A 



(28) 



then 



Ii^x) =-—— — = / dTj Xjij (x,T), 

OAji{X) Jq 

£«/+(*'))) 



{g(x,x)) 



lim 



(29) 



(30) 



This differs from the formula for (g) earlier by includ- 
ing a contact term, and which was avoided then by the 
condition on the sites i ^ j. 

With the choice (|27|) the WI (f25|) becomes 



D x «/(*)» = d x ((I(x)}) + i[A x , ((I(x)))} = 0. (31) 
The current I(x) has the same structure as A(x): 



I(x) = 



lW(x) I+(x) 



I-(x) I^(x) 
and the + component of (^Tj) is simply 

d x (/a + I^-I^A++A^I+-I+A^ 



(32) 



= 0. 



Differentiating with respect to A + 1 (x') and using (30) 
we obtain 

d x (g(x', x)) + i8(x - x') lim ((/^(z) - /&>(*))) = 0. 



A^0 



In the next section we will show that the second term 
here vanishes in all orders of perturbation theory. As a 



consequence, the mean conductance is independent of the 
positions of cross section x. Similarly, it is independent 
of x' as well: 

d x (g(x,x'))=d x ,(g(x,x'))=0. (33) 

This can be used to average (g) over these positions: 



1 



dx / dx'(g{x,x')) 



(34) 



lim 

L 1 n^O 





L pL 

dx / dx' 
o Jo 



5A+(x) 



+ I n (x)I+(x') 



Similarly, the second moment is obtained by taking 
four derivatives: 



(g(x 1 ,x' 1 )g(x 2 ,x' 2 )) 



lim 



<5 4 Z[A] 



IZ" SA+ 1 (x 1 )SA n (x[)SA+ 2 (x 2 )6A 22 (x' 2 ] 



(35) 



Again we can show that this expression does not depend 
on any Xi, and we write 

1 f^ 1 f f f 

var g — — j lim / dx\ I dx[ / dx 2 I dx' 2 

L A^o ^° J ° Ja J ° 

51+ (x[) 5I+(x> 2 ) | SI+(x' 2 ) 6I 22 (x 2 ) 

SAf^Xi) 8A£ 2 (x 2 ) SAf^xx) 8 A^fa) 
5I 22 (x 2 ) 51+ 2 {x' 2 ) 
' SA+ix!) 8A^(x[) 

-^I, MX2 ) + (1 



2) 
2) 



-^plX l{x [)I^ 2 ) + (+ 
& A n ( X V 

- In ( x i ( x 'i ) 7 22 t x i)I 22 (4 ) 



(36) 



where the subtraction of (g) 2 leaves only connected terms 
in the above expression. 

Now we specify the gauge potential even further, leav- 
ing only the A+ components nonzero, which is all we need 
to calculate moments of conductance. With this gauge 
choice we make the substitution ( ^o|) in the gauge invari- 
ant action S[A]. The result is S[A] = S + S int + S^iAJ, 
where So and Si nt are the same as before, Eqs. ([l7],[l8]), 
and 

S'[A X ] = -A Ps dx drtr (iA~ {d x tf - 26 t d a! && t ) 
Jo Jo 

+ iA+ (d x (b + btfb) - 2bd x tfb + 2btfd x btfb) 
+ A+A- + A+A-tfb + A+btfA- 
+ 2A + btfA-tfb - A-tfA-tf - A+bA+b 
-2A+bA+btfb - A+btfbA+btfb) . (37) 
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The currents I x and their derivatives necessary for the 
calculation of (g) and var g follow from ( p9| ) and (^) (we 
assume notation r = (x,t), r' = (x',t'), etc.): 



V x )Ia=o 
8A H (x) 



8A+(x) 



4i Ps / dr(d x tf -2tfd x btf) (r), 



= -4ip s / dr + &6+6) - 26^ c 6 t & 
Jo 

,c 

= -4p s 8(x - x') I dr (8 ik 8ji + (b^b) ik 8ji 



+8 lk {bb% +2{tfb) lk {bb%) (y), 
= -4p s S(x - x') I dr (8 ik 8ji + (btf) ik 5ji 



+8 ik (tfb) lj +2(btf) ik (tfb) lj ) (r), 

= 8p s 8(x - x') f dr (b* ki b* jt ) (r), 
Jo 

= 8p s S(x — x') I dr (b ik bij + b ik (bb'b)ij 
Jo 

+ (bb^b) lk b lj + (66t6) 4fc (66t6), J ) (r). (38) 



If we use all four blocks of A x in this derivation, wc 
would also obtain the following expressions for the diag- 
onal currents, necessary to prove Eq. (|33|): 



rC 

= -4ip s / dr (d x b% - tfd x b 

A=0 Jo 

-2tfd x bb%. (r), 



A=0 



Aip s I dr (d x bl 
Jo 

-2&&td x &&t) (r) 



(39) 



Substituting expressions (J3q) into Eqs. (34 39), we no- 



tice that the total derivative terms in give zero after 
integration over x's, thanks to the boundary conditions 
(|l9|). Dropping these terms, we obtain the exact expres- 
sion for the mean conductance 



4C Ps , ip 



(ff) = —ff + JJ I dx I dr 



o Jo 



«(& t &)ii + (» t )n+2(&t&) u (&&t) u » 

64p 2 II r L r c + + 

Hs 11 1 dx dT{bd x tfb-btfd x btfb) n (r) 



L 2 



)o Jo 



x dx' dT'(tfd x >btf) n (r') 
Jo Jo , 



(40) 



(from here on, the limits n — > and A — > will be 
implicit in expressions for observables such as (g) and 
varg). The expression for var g is more complicated and 



we do not reproduce it here. It will be clear from the 
next section, where we discuss the details of perturbative 
calculation, that only two terms in the above expression 
( pl0| ) contribute to the function <& of Eq. (0). Similar 
simplifications are even more drastic for the variance of 
the conductance. 



V. DETAILS OF PERTURBATIVE 
CALCULATION 

We have now assembled all the tools required to com- 
pute the mean conductance and its variance, and are 
ready to perform the explicit computation of the uni- 
versal functions <& and 

The free action So is diagonalized by Fourier trans- 
formation of the fields 6(r) to momentum space b(r) = 
C^ 1 Y. k , m h (v)4>k{x)e~ WmT , where p = [k,iu m ) and 
<f>k(x) = {2/L) 1 / 2 s\nkx: 



S ° = ^Y1 ( 4 ^ fc2 " 2iM o^m) tr&t(p)&( p ). 



(41) 



Periodic boundary conditions in the r-direction and the 
homogeneous ones in the x-direction ( |l9| ) imply that the 
frequencies oj m and momenta k take the values 



Lu m = 2irm/C, m — 0, ±1, . 
k d =nl/L, 1 = 1,2,... 



(42) 



The action (^lj) implies that the bare propagator of the 
bosons b (analogous to the diffuson of the usual isotropic 
metallic systems) is diagonal in replica indices: 



(( & pg(p)Mp')))n = <VV<W d 0(p), 



where 



do(p) = 



C 



CL 2 



X 



Ap s k 2 — 2iMouj m 47r 2 p s XI 2 — 2mm 



(43) 



(44) 



in which A = 2n 2 L 2 /L 2 as before. The subscript on 
the functional average indicates that it is taken using the 
action So (which does not contain A). The replica in- 
dex structure of the interaction vertex (|l8| ) is such that 
this diagonal property holds in all orders of perturbation 
theory, so that the full propagator is also diagonal when 
A = 0, ((b* q b rs ))A=o fx 5 pr 5 qs . Then it immediately fol- 
lows that the contributions to (g) quadratic in b disap- 
pear in the replica limit: (((&'' &)ii))a=o = ^Pa^n)) A=o oc 
8 a = n -> 0, and, similarly, (((fcfo t )ii)) A =o -> 0. 

The same argument applies to the first two terms in 
the expressions ( p^ ) for the diagonal currents. Moreover, 
it is easy to see that the quartic terms there after aver- 
aging also contain at least one factor of n, in all orders 
of perturbation theory. Then it follows that 
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lim 

n^O 



hm 

n^O 



0, 



(45) 



which proves the divergencelessness of the conductance, 
The bare propagator in real space 



d (r,V) 



L A sin fc;xsin k t x' iWm(T _ T <) 

2tt 2 C Ps ^ XI 2 - 2mm 6 



(46) 



has the dimensionful prefactor The interaction ver- 
tex S'int ( |l8| ) has the dimension Cp s /L oc £/L by inspec- 
tion. Then in the M-th order of perturbation theory the 
contribution to a general correlator of 2N bosonic oper- 
ators in real space, having 2M + N bare propagators in 
it, will be proportional to (L/£) M+N . Therefore to or- 
der LIE, which is all we need to obtain the function $ of 
Eq. (||), we have to calculate only two correlators in the 
zeroth order of perturbation theory. 
The first correlator is 



«(6 t 6)n(» t )ii» - to 



d (x,T + ;x,r) 



(47) 



(we used Wick's theorem and dropped the term disap- 
pearing in the replica limit), and the corresponding con- 
tribution to the function <f> is 



/ 



E 



XI 



2 — 2irim I 



\ e *> m3 n 



2 ^— ' XI 2 — 2m'm\ XI 2 — 2m'm2 



(48) 



where rj is a positive infinitesimal. When the frequency 
summations over m are done, we obtain the first two 
terms in Eq. (Q) . The second correlator we need is 



{{(bd x b^b) 11 (r)^d xl bb^ 11 (r'))) Q 
= d x d x ,d (r;r')d 2 (r';r), 



(49) 



and the corresponding contribution to $ gives the last 
term in Eq. (§ (using f(x)-f (y) = f(x)f(-y)/f(x-y)). 

The calculation of the variance of the conductance pro- 
ceeds along the same lines, except that the power count- 
ing and the replica limit leave us only with two correlators 
to compute. The first one is 

«(&t&)n(r)(66t) 22 (r') + (66 t )ii(r)(6t 6 ) 22 (r')» 

= 2d (r;r')do(r';r), (50) 

and the second is 

(((6 21 6 12 )(r)(6^6* 1 )(r'))) =^(r';r)- (51) 



The function $ is given then by 

„2 i-h i-C rL 



$ (L/L ) 



dx dr dx I dr' (2dg(r';r) 



L A Jo Jo Jo ./i) 
+d (r;r')d (r';r)) 



(52) 



An expression of this form is quite standard in the uni- 
versal conductance fluctuation theories, see for exam- 
ple, Eq. (5.3) in Ref. Q. The only difference is the 
anisotropic form of the diffusion propagator do(r; r'). In 
momentum space we obtain 

~ 2A 2 ^ / 2 1 

( 1 °>~ n 4 f^\\Xl 2 ~2mm\ 2 + {XI 2 -2mmf 

(53) 

Performing frequency summations we arrive at the ex- 
pression (||) [using f(x) + /(— x) = —1 and f'(x) — 
f(x)+f 2 (x) = -f(x)f(-x)i 

The limiting values (|9|-[T2|) are obtained using the 
asymptotic forms of the Bose function: 



l/x, x — > 0, 
e~ x , x — > +oo, 
— 1, x — > — oo. 



(54) 



In the ID limit L > Iq (A C 1) the expression 
reduces to 



1 °° 1 
9 E Ti 



I 4 



h,h = 



1 



1 Ifih + h) 2 



(55) 



These sums are easily evaluated, and we obtain Eq. (g). 
Similarly, the expression (|J) in the same limit L ^> Lo 
simplifies to 6tt~ 4 J2i^~ 4 = 1/15. In the opposite 2D 
chiral metal limit L <C L Q (A S> 1) all the terms in Eq. (Q) 
are exponentially small, and the largest of them gives 
Eq. (11). In contrast, the first term in Eq. (j^) diverges 
in this limit giving Eq. p2[). 

The two terms in eq. (|52|) correspond to two distinct 
correlation functions as above, or to two Feynman dia- 
grams. Mathur || evaluated only the first of these, in the 
two limits A — > oo and A — > 0. It happens to give the full 
contribution asymptotically as A — + oo. The remaining 
diagram is nonzero in general, and for A <C 1 contributes 
exactly 1/2 as much as the first one, to give the total in 
eq. A 



VI. CONCLUSION 

In conclusion, we considered the directed network 
model of edge states. In previous studies [@||@] it 
was shown that in a continuum limit this model can 
be mapped to a ID quantum ferromagnetic spin chain. 



9 



Three regimes, 2D chiral, ID metallic, and ID local- 
ized, separated by smooth crossovers, have been identi- 
fied for the model, and the universal crossover functions 
for the mean and variance of the conductance have been 
obtained for the crossover between ID regimes in our 
previous paper B. In this paper we use spin- wave per- 
turbation theory to obtain the corresponding universal 
functions for the other crossover, between the 2D chiral 
and ID metallic regimes. The results and their asymp- 
totic forms are given in Sec. H see Eqs. (@-[l|), and were 
discussed there. 
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